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Abstract 

S100A4 interacts with many binding partners upon Ca'^^ activation and is strongly associated with increased metastasis 
formation. In order to understand the role of the C-terminal random coil for the protein function we examined how small 
angle X-ray scattering of the wild-type S100A4 and its C-terminal deletion mutant (residues 1-88, A13) changes upon Ca^^ 
binding. We found that the scattering intensity of wild-type S100A4 changes substantially in the 0.15-0.25 q-range 
whereas a similar change is not visible in the C-terminus deleted mutant. Ensemble optimization SAXS modeling indicates 
that the entire C-terminus is extended when Ca^"^ is bound. Pulsed field gradient NIVIR measurements provide further 
support as the hydrodynamic radius in the wild-type protein increases upon Ca^^ binding while the radius of A13 mutant 
does not change. IVlolecular dynamics simulations provide a rational explanation of the structural transition: the positively 
charged C-terminal residues associate with the negatively charged residues of the Ca^Mree EF-hands and these interactions 
loosen up considerably upon Ca^^-binding. As a consequence the A13 mutant has increased Ca^^ affinity and is constantly 
loaded at Ca^^ concentration ranges typically present in cells. The activation of the entire C-terminal random coil may play a 
role in mediating interaction with selected partner proteins of S100A4. 
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introduction 

Through metastasis cancer cells migrate from the primary 
tumor to distant organs to establish secondary growth and this is 
considered the most common cause of death amongst cancer 
patients. The transformation of cancer cells to a metastatic 
phenotype involves a variety of extracellular and intracellular 
alterations such as cell signaling, interactions with the extracellular 
matrix (ECM) and cytoskeleton remodeling. S100A4 (also known 
as mtsl, metastasin, calvasculin, p9Ka) has been strongly 
associated with increased metastasis properties in several animal 
models [1,2]. 

The interaction of S100A4 to different binding partners (for 
example p53, annexin A2, CD 16, non-muscle myosin II) is 
mediated by Ca """-induced conformational changes. It has been 
widely recognized that the conformational change upon Ca^"*"- 
binding opens up a hydrophobic cleft on the protein surface [3-6] . 
Recentiy a series of X-ray and a NMR structure of S100A4 
complexed with a NMIIA fragment revealed detailed information 
about the binding sequence and the molecular interactions 
between S100A4 and the NMIIA fragment [7-9]. 



Mutants with shortened C-terminal regions [10,1 1] reduced the 
ability of S 1 00A4 to bind to NMIIA and to promote metastasis. 
The deletion of the last two C-terminal amino acids as well as the 
last fourteen amino acids of S 1 00 A4 reduced metastasis formation 
compared to the wild-type S100A4 [10]. Metastasis was 
undetectable after the deletion of the 3, 4 and 5 C-terminal 
amino acids of S 1 00A4 in the same in vivo study, while as little as 
the removal or mutation of the last two residues reduced the 
binding affmity to NMIIA by an order of magnitude [10]. The 
crystal and NMR structures determined so far did not provide any 
clue to rationalize this observation, as the conformation of the C- 
terminus does not appear to be well defined. As such the C- 
terminus of S 1 00A4 is one of the most conformationally variable 
regions in the crystal [3,8,12] structures. Through NMR studies of 
wild-type S100A4 [13,14] chemical shift differences at the stem of 
the C-terminal random coil upon Ca^"^-binding were detected, 
which was interpreted as an extension of helix 4 by approximately 
4 residues (residues 86-90) whereas for the last 1 1 residues a 
significant change in chemical shifts was not visible [13,14]. The 
relaxation parameters of Ca "^-free S100A4 indicated slow 
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dynamics in tliis region, but this has not been compared to Ga^*- 
bound or target bound states of S100A4. [15]. 

The long C-terminal random coil region is characteristic to the 
S100A4 and SIOOAIO proteins, the closely related SlOOAl, 
S100A6 and SIOOB proteins have shorter C-terminus. It has a 
hydrophobic stem (Phe-89, 90 and 93, Pro-94) and a tip populated 
by bulky basic residues (Lys-96, 100 and 101, Arg-99) (Figure 1). 
In two crystal structures of Ca^"^-bound peptide-free S100A4 the 
hydrophobic cleft of each subunit is occupied by the hydrophobic 
stem of the C-termiuus of an adjacent dimer [6,12], which may 
explain the formation of S100A4 oligomers that were observed 
extraceUularly [16]. This cleft is occupied by the NMIIA peptide 
in the S100A4:NMIIA complex [7,8] and the hydrophobic stems 
interact with one other in the crystal forming an elaborate 
hydrophobic knot between two dimers [8]. Additionally, this 
unusually long C-terminus with the positively charged tip has a 
crucial role in annexin A2 (ANXA2) mediated tissue plasminogen 
activator activation which leads to plasmin formation and finally 
MMP activation. [5,17]. 

Since the role of the C-terminal random coil is not well 
understood in the mechanism of Ca^^- binding we investigated 
wild-type and truncated S100A4 (A 13 mutant truncated to residue 
88) by solution small-angle X-ray scattering and NMR spectros- 
copy and rationalized the observed changes by molecular 
dynamics simulations. We also examined how the deletion of the 
C-terminal region affects Ca^"'"-binding by isothermal titration 
calorimetry and solved the crystal structure of A13 mutant in 
complex with the NMIIA interacting peptide. 

Materials and Methods 

Production and Purification of the S100A4 Mutants and 
the Myosin IIA Fragment 

The S100A4 variants and NMIIA fragment were obtained as 
described in our recently published work [8]. Briefly, the His- 
tagged human S100A4 variants were expressed in E. coli BL21 
strain, and purified by Ni^"''-afrmity chromatography. After 
cleavage of the Hise-tag with Tobacco Etch Virus protease, the 
sample were applied to phenyl-Sepharose column, washed 
thoroughly and eluted with EGTA containing buffer. The purified 
S100A4 proteins were dialysed against 20 mM Hepes pH 7.5, 
20 mM NaCl, 0.1 mM TCEP buffer, concentrated and pooled at 
— 70°C. Prior to SAXS measurements, the S100A4 alone or in 
complex with MPT peptide were apphed to Superdex 75 10/300 
GL gel filtration column (GE Healthcare, Little Chalfont, UK) 
and the peak corresponding to ~24 kDa molecular weight was 
collected. The NMIIA fragment MPT (Tyr-Argl894-Lysl937) 
was expressed with an N-terminal Hisio-ubiquitiii fusion in E. coli 

Protein Sequence 



80 90 100 

S100A4 cvFLsciAMMCNEiaaEGiagpiaop iaaa loi 

SlOOAl VVLVAALTVACNNFFWENS 94 

S100A6 VTFLGALALI YNEALKG 90 

SIOOB MAFVAMVTTACHEFFEHE 92 

SIOOAIO FSLIAGLTIACNDYFVVHMKQKGKK . 97 



^hydrophobic stem 
^basic residues 

Figure 1. Sequence alignment of several S100 protein family 
members. The C-terminal region of SlOOAl, S100A6 and SIOOB is 
compared to the last 25 residues of S100A4. The C-terminus of S100A4 
consists of a hydrophobic stem and basic tip residues. 
doi:1 0.1 371/journal.pone.0097654.g001 
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BL21 strain, and purified by Ni'^-affinity chromatography. The 
His-tagged ubiquitin was removed by Yeast Ubiquitin Hydrolase, 
while it was dialysed against buffer containing 20 mM TRIS (pH 
8.0), 100 mM NaCl and 0.2 mM DTT. The completely digested 
sample was applied to Ni^''"-affiiiity column. The peptide was in 
the flow-through, which was fmaUy purified by reverse-phase 
HPLC on a Jupiter 300 C18 column (Phenomenex, Torrance, 
CA). 

Isothermal Titration Calorimetry 

Titrations were carried out in 20 mM Hepes pH 7.5, 150 mM 
NaCl, and 1 mM TCEP using a Microcal VP-ITC apparatus (GE 
Healthcare, Little Chalfont, UK). 100 ^M S100A4 variants were 
titrated with 3 mM CaClj at 310 K. This temperature is optimal 
to resolve the two binding sites as described previously [6]. 39 
injections were performed with 350 sec time intervals between 
injections. The Origin for ITC 5.0 (OriginLab) software package 
was used for data processing and the model "Two Set of Sites" 
was fitted. Peptide binding experiments were carried out at 298 K 
as described previously [8]. 

SAXS Measurement and Data Analysis 

SAXS data were collected on the beamlines ID 14-3 and BM-29 
at the European Synchrotron Radiation Facilities (ESRF, 
Grenoble, France). Data from protein solutions at a concentration 
range from 5 mg/ ml to ~ 1 mg/ ml were collected under 
continuous flow with an X-ray wavelength of 0.931 A at 283 K 
and an exposure time of 7 s per frame at ID14-3, and 0.992 A and 
an exposure time of 1 s per frame at BM-29. 10 frames for each 
concentration were recorded (Detector: PUatus 1 M (Dectris), with 
0.01-5 nm ' range). Scattering curves measured at different 
concentrations were merged and processed in Primus [18]. 
Theoretical scattering profiles of the high resolution structures 
were calculated and compared to experimental scattering profiles 
using the program Crysol [19] with default parameters. Ensembles 
of S100A4 with constant core domain (residues 1-85) with variable 
C-terminal region were generated using the program EOM [20]. 
The calculated intensity from the subset of these ensembles was 
optimized against the Ca^^-free and bound experimental data. 

NMR Measurements and Data Analysis 

NMR spectra were collected on a Bruker Avance III 
spectrometer operating at 700.17 MHz frequency for 'H, 
equipped with a 5-mm 'H/''^C/''^N probe-head with z-gradient. 
Typical sample composition was 1 mM protein, 150 mM NaCl, 
20 mM Hepes, 2 mM TCEP, 5 mM EGTA, 10% DjO. Data 
were collected on the sample containing the Ca^^-free protein, 
followed by addition of Ca^"^ solution and the resulting Ca^"*" 
concentration was 9 mM. Measurements were done at 287 K, and 
the temperature was calibrated with standard methanol solution. 
Diffusion measurements were done with a standard Bruker pulse 
sequence using stimulated echo, bipolar pulses and 3-9-19 
sequence for water suppression. The gradient calibration constant 
was determined using doped H2O sample at 298 K. The sample 
temperature was stable within 0.01 K during one experiment. ID 
'H spectra were recorded in pseudo 2D mode and gradient 
strength was varied linearly in 32 steps between 2% and 95% of its 
maximum value. Delays were optimized in order to achieve a total 
decay in protein signal, thus the apphed diffusion time (A) was 170 
and 200 ms, and gradient length (5) was 4 and 5 ms. For each 
sample at least 3 measurements in 8 K data points with 128, 256 
or 512 transients were acquired. Data processing and analysis was 
done in TopSpin. Upon processing only zero order phase 
correction was applied and integration was performed in regions 
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Figure 2. Scattering curves of Ca^Mree, Ca^*-bound and iVIPT-bound S100A4. Differences in tine WT [cyan) and the Ca'^^-bound WT {red) 
are mainly observed at 0.15 A^'-0.25 A^' (A). The scattering profile of the A13 mutant {purple) changes less upon Ca^*-binding {orange) (B). Typical 
one-dimensional NMR spectra for the studied systems acquired at 700.17 MHz, zoomed to the aliphatic proton region. From bottom to top: buffer; 
WT; WT-Ca^*; A1 3; A1 3-Ca^* (C). Typical examples of EOM [20] extended (D) and compact SI 00A4 models (E). Distribution of the radii of gyration of 
the generated model ensemble {black) for ensemble optimization method [20], the best ensemble fitting the Ca^*-free {cyan) and Ca^*-bound 
{magenta) WT S100A4 SAXS data (F). Scattering curve differences between the MPT-bound WT {red) and the MPT-bound A13 mutant {orange) are 
more diffuse (G). 

doi:1 0.1 371 /journal.pone.0097654.g002 

where no buffer peal^s were present. In this respect two chemical 
shift ranges were chosen: (1.33-1.15) ppm and (0.82-0.20) ppm. 

The decay of the protein signal integrated intensity with 
gradient strength was fitted to a single component. This approach 
was satisfactory in all cases - meaning homodimers are present in 
the solution. However, in the case of the Ca^"*"-bound WT protein 
sample the possible presence of oligomers was checked by fitting 
two components. Even in this case one component, corresponding 
to a homodimer, was sufficient to model the data. Diffusion 
constant calculation was based on fitting the variation of integrated 




0 10 20 30 40 50 60 70 80 90 100 

Time (ns) 

Figure 3. Radius of gyration of 100 ns molecular dynamics trajectories. Ca^*-free S100A4 WT {black) and AAA constructs {red) are compared 
with snapshots of the structure at 10, 30, 50, 70 and 90 ns the C-terminal random coil, the N- and C-terminal EF-hand are colored yellow, magenta and 
cyan, respectively (A). Trajectories of the same constructs when in the Ca'^^-bound form (B). 
doi:1 0.1 371/journal.pone.0097654.g003 



intensity values I as function of gradient strength G; according to 
the equation: 

where D is the diffusion constant, A is a constant containing 
diffusion time, diffusion gradient length values and the calibrated 
gradient constant. Fitting was done on 29-31 points. From the 
obtained diffusion constant values, assuming aU molecules are 
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Figure 4. Thermodynamics of Ca^'^-binding to wild-type (A) and A13 (B) S100A4 measured by isotKiermai titration calorimetry. Plots 
represent the raw heat effects as a function of time and the binding isotherms as a function of the molar ratio of Ca^* to the protein (100 \iM). 
Calculated thermodynamic parameters are shown in Table 2. 
doi:1 0.1 371/journal.pone.0097654.g004 



spherical, hydrodynamic radii were calculated on the basis of the 
Stokes-Einstein equation: 



updated every 2 ps and each simulation was carried out for 
100 ns. 



Rk- 



kT 
6nr\D 



where k is the Boltzmann constant, T the absolute temperature 
and ri is the viscosity of the aqueous solution at the given 
temperature. For each measurement two diffusion constant values 
were calculated for the above mentioned two H chemical shift 
regions, thus the given values are the average of at least 6 data 
points. A typical decay curve is shown on Figure SI and the 
determined diffusion constants and corresponding Rj, values are 
listed in Table SI. 

Molecular Dynamics Simulations 

A standard GROMACS protocol was used for construction of 
the systems in Table S2. Protein molecules were put in a 
dodecahedron simulation box and solvated in SPC/E water. For 
all systems, the necessary number of counter ions was added to 
neutralize system charge and additional sodium and chloride ions 
were provided until approximately 0.1 M concentration was 
achieved. The OPLS-AA/L force field was used for all simula- 
tions. All systems were energy minimized and subsequently heated 
up to 300 K temperature. In all simulations a time step of 2 fs was 
used. The simulations were performed in the isobaric-isothermal 
ensemble (NPT) and the temperature and pressure (1 bar) was 
scaled by the modified Berendsen thermostat [21] and the 
ParrineUo-Rahman barostat [22] with 0.1 ps and 2 ps sampling, 
respectively. The van der Waals and electrostatic interactions were 
truncated using the Verlet cutoff scheme [23]. The trajectory was 



Results and Discussion 

Small angle X-ray scattering (SAXS) and pulsed field gradient 
NMR spectroscopy were applied in order to study the structural 
effect of Ca^"'"-binding and the hydrodynamic radii (Rh) for the 
S100A4 WT and A13 mutant (see Table 1 for notations) in their 
Ca^'*"-free and Ca^^-bound state (Figure 2A, 2B and 2C). While the 
scattering intensity of the Al3 mutant did not change substantially 
upon Ca^-'-binding (Figure 2B), the wUd-type S100A4 SAXS 
intensity increases noticeably in the 0.15-0.25 A ' q range. In the 
protein concentration range (35-113 |iM) used in this SAXS 
experiment the Guinier plot [24] is linear (Figure S2) and the 
experimental Rg values obtained at the lowest concentrations are 
summarized in Table SI. Modeling an ensemble containing 
random conformers of the C-terminal region [20] yielded distinct 
shapes such as shown on Figure 2D and 2E. In the Ca^'''-bound 
form an extremely extended conformation of the C-terminal 
region dominates the modeled ensemble (Figure 2D and 2F) 
(average Rg 22 A, maximum size 94 A, 0.81) whereas in the 
Ca^"'^-free form (Figure 2F) there is a significant peak at the 
approximate Rg of 1 9 A representing a more compact conforma- 
tion (Figure 2E) (average ensemble Rg 21 A, average maximum 
size 85 A, 0.89). A similar tendency is shown in the 
hydrodynamic radii (RiJ between the Ca "^-free (25.6±0.4 A) 
and bound form (33.1 ±1.6 A) of S100A4, which were calculated 
from the diffusion coefficients obtained by pulse-field gradient H- 
NMR experiments. On the other hand the Rj, of Al3 mutant does 
not change substantially from Ca^''"-free (25.8±0.3 A) to Ca^"""- 
bound form (24.8±0.6 A). The aliphatic region of the 'H-NMR 
spectra of the four systems is shown on Figure 2C. The 'H spectra 
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indicate broader peaks for the WT protein compared to the A13 
mutant. While tlie addition of Ca^"*" ions to the Al3 mutant does 
not alter the signal intensities, an intensity decrease is observed for 
the WT protein. This is in line with the formation of oligomers 
that possess peaks broadening below the detection limit. Data 
analysis of the diffusion measurements was performed in regions 
where only protein signals are present (Figure 2C). For the Ca^^- 
free WT and A 1 3 mutants the Rh values do not differ significantly, 
indicating both proteins have similar size. This means the longer 
WT molecule has to be more compact than the shortened mutant. 
The addition of Ca^"*" ions causes a small decrease in the Rh value 
for the A13 mutant, whereas for the WT protein a considerable 
increase in size is obtained in the presence of Ca^""" in line with a 
putative extended conformation. 

The relationship between the obtained hydrodynamic radius 
and radius of gyration for spherical molecules is Rg= (3/5)'''^ Rj, 
[25]. The ratio of 0.77 agrees remarkably well for the Al3 mutants 
with R^ estimated to be 19.4±0.8 Aand 19.1±1.0 A for the Ca^"^- 

2+ 

free and Ca -bound form, respectively based on the Guinier 
equation. Even a corresponding compression of the R^ can be 
observed upon Ca^"*^ binding (Rh of 25.8±0.3 A and 24.8±0.6 A, 
respectively). For the WT protein the agreement is less perfect, 
indicating that the shape is different. However, the Rg of the 
smaller component from the optimized ensemble of the Ca^"'"-free 
model (Rg approximately 1 9 A) is in good agreement with the Rj, 
of 25.6±0.4 A and .similar tendencies in size changes obtained 
from two independent methods illustrating a specific structural 
role of the C-terminal random coU. 

The agreement of the high-resolution structures to the scattering 
intensities was also evaluated using the program Crysol [19] as 
described in the Materials and Methods. 

The values of the fits of Ca^"'"-free WT solution scattering 
compared to the Ca^'^'-free WT NMR structure model ensemble 
(PDB code 1M31) are in die range of 1.85-1.56 (Figure S3A and 
S3B). For the Ca^'^-free A13 mutant the %' values are between 
1.19 and 1.54, but artificial truncation of the model to compose 
only residues 1-88, on the other hand, provided a much better fit 
against the A13 mutant data with a range of 0.90-1.09 (Figure 
S3A and S3C). The Ca^'^-bound S100A4 SAXS intensity 
prediction from the crystal structures 3C1V [12], 3CGA [3] and 



2Q91 [6] PDB entries are comparable, showing an agreement of a 
1.05, 1.19 and 1.04, respectively (Figure S3D). 

The SAXS intensity curve of the A 1 3 mutant is very similar to 
the Ca^"''-bound and free form, which may suggest that either their 
structure is similar or the limited resolution of the SAXS technique 
is unable to distinguish clearly these functionally distinct states. On 
the other hand, if the structures are similar, the A 1 3 mutant may 
be expected to bind MPT with similar affmity in the Ca^"'"-free 
state as in the Ca^"'^-bound state. Figure S4 and Table S3 show that 
Ca^'*"-bound Al3 S100A4, just as the WT protein displays similar 
affmity to MPT close to the detection limit of the ITC technique, 
while the Ca2"^-free state of A13 S100A4 does not bind MPT m 
accordance with the previous findings with WT S100A4 [26]. This 
finding indicates that the more pronounced SAXS intensity 
differences between the Ca^'*"-free and bound states of WT 
S100A4 is related to the differences in the conformational 
ensemble of the flexible C-terminal random coil. Since the ITC 
measurements did not reveal the true affinity of MPT to S 1 00A4, 
we also performed fluorescence polarization measurements. These 
showed that WT and A13 S100A4 have similar affmity to 
fluorescein labelled MPT with subnanomolar of 
0.12±0.06nM (SEM) and 0.28±0.11 nM (SEM), respectively 
(FigTjre S5A, B). Similarly, the disassembly of dimeric NMIIA rod 
filament (fragment 1712Q;1960E) by S100A4 [8] is not afffected 
substantially by the C-terminal truncation (Figxire S5C, D). 

When a peptide derived from NMIIA (MPT: residues Argl894- 
Lysl937) was added to the Ca^'^^-activated WT and the Al3 
mutant the scattering intensity of both variants change similarly 
(Figure 2G), differences between the wild-type and A13 complex 
are more spread out and difficult to interpret without additional 
information. Therefore we solved the crystal structure of Ca^"""- 
saturated A13 mutant in complex with MPT peptide (Supporting 
information. Figure S6, Table S4), which apart from some 
differences observed (Supporting information. Figure S7), essen- 
tially confirmed the asymmetric binding mode of MPT observed 
in complexes of the full length protein [7,8] and A8 C-terminal 
truncated mutant [9]. 

To shine light on the mechanism of Ca^''"-activation in the C- 
terminal random coil molecular dynamics simulations were 
performed with the GROMACS package [27] on six starting 



Table 1. Description of the proteins and complexes 


compared in the text 


and figures. 






Notation 


Experimental data 


Substitutions Deletion 


High resolution structures 


Ca^+- and MPT-bound F45WSer 


F45WSer 


X-ray 


C3S, C81S, C86S, F45W 


Ca^+- and MPT-bound A13 


A13Ser 


X-ray 


C3S, C81S, C86S F89-K101 


Ca^+-bound WT 


3C1V 


X-ray[12] 




Ca^*-bound A8 


4HSZ 


X-ray[9] 


F93-K101 


Ca^+- and MPT-bound WT 


2LNK 


NMR[7] 


- (Ca^^ ions not modeled) 


Ca^*-free WT 


1M31 


NMR[14] 




Solution scattering and NMR data 


WT 


WT 


SAXS/NMR 




Ca^+-bound WT 


Ca^+-WT 


SAXS/NMR 




A13 


A13 


SAXS/NMR 


F89-K101 


Ca^+-bound A13 


Ca^+-A13 


SAXS/NMR 


F89-K101 


Ca^+-and MPT-bound WT 


Ca=^+-WT 4-MPT 


SAXS 




Ca^+-and MPT-bound A13 


Ca=^+-A13-hMPT 


SAXS 


F89-K101 


doi:l 0.1 371 /journal.pone.0097654.t001 
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configurations as presented in Table S2. When Ca ^ ions were not 
present the C-termini rapidly became locked on the surface of the 
core molecule, predominantly associating with the negatively 
charged residues of the EF hands (Movie SI). As Figure 3 shows 
similar stabilization was not observed in the Ca^"*" bound form. 
Here the C-terminus remained dynamic even after 50 ns. This 
suggests an electrostatic control of the dynamics of the C-terminal 
random coU: when Ca^"*" binds to the EF-hands the Coulombic 
interaction between the EF-hands and the C-terminus weakens 
and the C-terminus gains conformational freedom (Movie S2). 
The trajectories may be biased by the starting configuration; 
therefore we repeated the simulation with the exact same starting 
coordinates as before except that the last three residues were 
replaced by alanines. In this case a curious reversal of the 
dynamical behavior of the C-terminus occurred: as Figxire 3 
shows, the Ca^''"-free trajectory indicates a more dynamic C- 
terminus whereas in the Ca^'*'-bound form the C-terminus 
becomes strongly associated by the hydrophobic groove of the 
Ca^'*"-bound, open form of S100A4. In the absence of positive 
charges the association with the empty EF hands is broken, but the 
alanines have a hydrophobic character and a sufficiendy small 
volume to bind into the open S100A4 groove which radically 
reduces the radius of gyration of the complex and renders it 
dynamically inactive. This suggests an additional role of the 
terminal residues: they prevent the autoinhibition by the C- 
terminus by having a hydrophilic character and their large side 
chain stericaUy interferes with the binding. The A 1 3 mutant both 
in the free and Ca^''"-bound form remained stable for 100 ns in 
molecular dynamics simulations. 

The calculated entropy difference using quasi harmonic 
approximation between the Ca^'''-bound and Ca^''"-free form of 
S100A4 is -0.60, - 1 .98 and - 1 .64 kJ/mol/K, for the WT, AAA 
and A13 constructs, respectively. Even though the diversity of 
conformations in the immediate vicinity of the Ca*""" ions is 
decreased as compared to the Ca^'''-free form, and in this respect 
the C-terminal variants do not differ, in the wild-type C-terminal 
region the conformational freedom increases upon Ca^'^'-binding 
leading to a reduced entropy difference between the two states. 

If the S100A4 C-terminus interacts specifically with the EF- 
hands, it should influence the Ca^'^-binding properties of S100A4. 
This is indeed the case: in the A13 mutant the Ca^"''-binding of the 
high-affmity site is exothermic (in contrast to endothermic Ca^""" 
binding in the WT protein) and the entropic component of the 
reaction is about one fourth of the value obtained for the wild-type 
protein (Figure 4A, B) (Table 2). The Ca^'^-binding property of the 
wild-type S100A4 is similar to previous observations [6]. This 
relatively high entropic component of Ca^'^-brnding could be 
explained by the detachment of the disordered C-terminus. 
Importantly, the Ca^''"-binding affinity of the high affinity site is 
forty times higher in the A 1 3 mutant than the affinity determined 
for the wild-type S100A4. Since the solvent effectively screens long 
range electrostatic interactions, the removal of positive charges 
upon the C-terminal truncation alone does not explain the 
increased affinity. The unstructured C-terminus of S100A4 has to 
specifically interact with the high-affinity EF hand of S 1 00A4 to 
reduce the Ca^"'"-aflBiiity and interfere with Ca'^-binding. The last 
1 1 residues of C-terminal random coil did not have a specific 
conformation in the Ca^''"-free state and switching between EF- 
hands occurred in the 1 00 ns simulation, illustrating why it may be 
difficult to use NMR chemical shifts to model Ca^'^-dependent 
local structural changes this region. [13,14]. 

With this high affinity, the C-terminally truncated S100A4 
would remain coiistitutively active inside of the cells and this alone 
would impair its Ca^^'^-dependent regulatory functions. This may 
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be of more broad ranging significance than putative binding 
affinity reductions [10] to individual interaction partners such as 
NMIIA and explain the strong influence of C-terminal truncation 
of SI00A4 in in vivo models. This is especially relevant since 
simplified models of NMIIA such as the shorter NMIIA derived 
peptide (1908-1923) [9], the MPT peptide in this study did not 
show significant reduction in binding affinity to truncated forms of 
S100A4 (A8 and A13 mutants, respectively). The A13 mutant also 
disassembled filaments of NMIIA fragment 1712Q,-1960E with 
similar efficiency as the wild-type protein. Nevertheless the 
dynamic Ca^"'^-activated elongated random-coil at C-terminus of 
S100A4 may be important for interacting with native NMIIA 
filaments. Binding to other partner proteins such as annexin A2 or 
p53 may also depend on an intact activated C-terminus in 
promoting metastasis in tumor cells, which require further studies 
[5]. 

The present results provide an insight to the dynamic 
mechanism of C-terminal random coil in S100A4 as a mediator 
of S100A4-driven metastasis and shine Ught on its role in tuning 
the Ca^"'"-binding affinity of S100A4. These results also suggest 
that locking the C-terminus to the core domain may be an 
alternative strategy of inhibiting the metastasis-promoting activities 
of S100A4. 

Supporting Information 

Figure SI Typical protein signal decay curve as 
function of the applied gradient strength. Points indicate 
measured values, while the continuous line is the fitted curve. 
(TIF) 

Figure S2 Guinier plots of the datasets that were 
merged for further data analysis and modeling. The 

Guinier region was extrapolated to the beam stop (red line). The 
data sets are plotted in rings with different colors representing 
different protein concentrations. 
(TIF) 

Figure S3 Comparison of the theoretical SAXS intensi- 
ties of the high-resolution structures with the experi- 
mental SAXS scattering. The experimental scattering curves 
are shown in black circles. The theoretical scattering curves were 
calculated with the program Crysol (21) and are shown in colored 
lines. The x^-values of the individual models (model number 1-20) 
of the Ca^Vree NMR structure ensemble (PDB code 1M31) [7] 
were compared to the Ca^^-free SAXS scattering curves (A). The 
calculated scattering cur\'e of the lowest energy NMR model 
1M31 against the Ca^'^-WT SAXS (B) and the Ca^'"-free A13 data 
(C), respectively. Truncated NMR models are indicated with T. 
Three Ca^"'-bound crystal structures (PDB code 3C1V, 2Q91, 
3CGA) are also compared to the Ca^"*^-bound WT SAXS 
scattering curve (D). 
(TIF) 

Figure S4 Thermodynamic analysis of S100A4A13- 
MPT interaction. 75 |XM S100A4A13 titrated with MPT at 
25°C in the presence of 1 mM CaCla (A) or 1 mM EGTA (B). In 
the absence of Ca^"*" no interaction was detected. Calculated 
thermodynamic parameters are shown in Table S3. 
(TIF) 

Figure S5 Fluorescence polarization measurements 
using (A) WT S100A4 and (B) A13 S100A4 and fluorescein 



labelled MPT peptide. Optical density changes at 320 nm 
when NMIIA 1712Q,-1960E rod fragment is titrated with (C) WT 
and(D) A13 S100A4. 
(TIF) 

Figure S6 Three-dimensional structure of the Ca^*- 
activated, MPT-bound A13Ser (A) and F45WSer (B) 
S100A4. Subunit A is shown in green, subunit B in blue and the 
bound MPT in yellow. Helices (H) and the N- and C-terminus (N, 
C) of the bound peptide are indicated. 
(TIF) 

Figure S7 Distance differences in the MPT-bound 
F45WSer and AlSSer. The superposition of the individucJ 

chains was performed with LSQIv'tAN [8] . Distance plot of subunit 
A and B of S100A4 (A) and the bound MPT peptide (B). The 
superposition of the three dimensional structure of chain A and 
chain B are shown in (C) and (D), of the bound peptide is shown in 
(E). 
(TIF) 

Table SI Experimental Rg values and hydrodynamic 
radius obtained by SAXS and NMR respectively. 

(DOCX) 

Table S2 Specification of the constructs used in differ- 
ent MD simulations. 

(DOCX) 

Table S3 Thermodynamic parameters of peptide bind- 
ing of the S100A4 variants determined by ITC measure- 
ments. 

(DOCX) 

Table S4 Crystallographic table. 

(DOCX) 

Movie SI MD trajectories of the Ca^^-free S100A4. The 

time interval visualized is between 20 ns and 100 ns. The last 17 
residues of the C-terminal region are shown in red and the EF- 
hand motifs in cyan. 
(WMV) 

Movie S2 MD trajectories of the Ca^^-bound S100A4. 

The time interval visualized is between 20 ns and 100 ns. The last 
1 7 residues of the C-terminal region are shown in red and the EF- 
hand motifs in cyan. 
(WMV) 

Text SI Supporting Online Information and Supple- 
mentary Methods. 

(DOCX) 
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